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time- dependent singular oscillator problem is solved in this approach. The important role of repre- 
sentation group theory is discussed in this context. 
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The method of disentangling expressions, containing non-commuting operators (FDM), suggested by Feynman in 
[l[ , gave an elegant solution of the harmonic oscillator excitation under the arbitrary time-dependent external force. 
In further developments the FDM was applied to some other non-stationary quantum mechanical problems, see 0, Q 
and [J] for review. And it was shown that the transition matrix elements calculation became much simpler if the 
FDM was supplied by some considerations from the representation theory of SU(2) or SU(1, 1) groups. 

In this paper we apply the FDM to a non-quadratic system — singular oscillator with a variable frequency w(t). 
We obtain the self-contained analytic expressions for the transition amplitudes between states with definite quantum 
i numbers (at t — > ±oc) and calculate the generating functions for transition probabilities. 

The problems considered in P, [|[ can be generalized to a model of a singular oscillator with variable frequency: 

£5 ■ < x < +oo, g — const, g > — 1, h = m = 1. 

a ; 

The frequency u)(t) is an arbitrary real time function. As usual, we propose the boundary conditions: 



O 

oo 



(t) — > oj± at t — > ±oo 



which allows one to define the final and initial states of the oscillator. 

It is well known that at a fixed t the instantaneous spectrum of the Hamiltonian ([T]) is equidistant (see, e.g., 0) 

: E n = 2u(n-j), J = -\-\Vl + ~9, n = 0,l,2,.... (2) 



We note further that the operators 



Ji = ife» + «p), J 2 = i(z 2 - P 2 )--^, j = 1 ^+ p ^ + ^ (3) 

4 4 v ' 16x z 4 v ' lbx z 

' satisfy the standard commutation relations of the su(l, 1) algebra: 

c3 ■ [^1,^2] = -iJo, [J2, M = iJi, [Jo, Ji] = iJz, (4) 

and the Hamitonian Q is a linear combination of operators Jo and J2' 

H{t) = (u; 2 {t) + 1) J + (u; 2 {t) - l) J 2 . (5) 

The instantaneous eigenfunctions of the Hamiltonian (JTJ) realize the irreducible unitary infinite-dimensional repre- 
sentation of the non-compact su(l,l) algebra. The corresponding Casimir operator ("angular momentum" squared) 
proves to be a constant and can be calculated directly: 



3> = J*-Jl-4 = ^=j{j + l) (6) 



so the weight of this representation is j. 
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For the simplest case w+ = lo- = 1 the initial and the final states are the eigenfunctions of the Jo operator: 

J ip„ = A„V>™, A„ = n-j. 

According to Ref. [3|], the transition amplitude between initial m) and final \n) states can be expressed in terms of 
the generalized Wigner function for the irreducible representation of the su(l, 1) algebra with weight j: 



Aj) 
J n—i 



n,m = 0,1,2,.... (7) 



The latter can be obtained by an analytic continuation of a standard Wigner function for the compact su(2) algebra; 
the details of this technique were described in the Appendix A in [3J] . In group theory such a method is known as the 
"Weyl unitary trick" . 

A generalization of this approach to the case of unequal initial and final frequencies oj + ^ lu- is quite obvious, as, 
according to ([5]), the transformation from the initial state basis to the final state one is simply a unitary rotation 
around the axis 1, i.e., an element of the quasi-unitary group SU(1, 1). Omitting intermediate calculations (all details 
can be found in [6(), we derive the following expression for the transition probabilities between states \m,u>-) and 
\n,uj+) (see Eq. (13) from Ref. for comparison): 

Wmn = {a -5)!)2g!r(g-2j-) pI " S(1 ~ P) ~ 2j [2 ^ l( "^ L ~ 2i; L ~ S + 15 P)]2 ' (8) 
where j is defined in © and 

L = max(m, n), S = min(rn, n), L — S=\m — n\. 

The p parameter, < p < 1, has the same sense as in the well-known problem of transitions in a regular oscillator 
(see, e.g., [7|, l8j) and can be calculated from the classical oscillator equation of motion. The formula ([8|) furnishes the 
ultimate answer to the problem given. 

The Gauss hypergeometric function 2-F1 in ((5J) has its first argument being integer and negative (or zero), so it 
reduces to the Jacobi polynomial. Making necessary transformations, we obtain from |8[): 



w 



run 



W. 



mn 



n! 1 (m — 2j) L 

- ^^r P m - n {l P)- 23 [p0»- ».-«-D(i - 2p) 2 
ml 1 (n — 2j) L 



n > to, 

(9) 

m > n, 



which coincides with the standard quantum-mechanical result for the transitions in the time-dependent singular 
oscillator obtained by means of the Shroedinger equation solution (see [9j and references therein). Consider some 
special cases of (JS]),©: 

a) The excitation of the ground (vacuum) state: to = 0, 



In particular, for g = — 1 we have 



--2fro'' , < 1 -'>-"■ (10) 



■i, m*> = p"(i-p). (ll) 



This is the utmost value of the g constant, as at g < — 1 the particle is collapsing into the potential centre [5[, so the 
Hamiltonian is not the self-conjugated operator any longer. It is significant that the expression (fTT|) corresponds to 
the distribution of charged scalar bos on p airs, created in vacuum with a definite momentum p under the influence of 
a variable homogeneous electric field jldl. fll|. 

b) At g = (i.e. j — —3/4) the problem considered corresponds to the regular oscillator with time-dependent 
frequency on the positive half-axis; in fact, the formulas © for this case should coincide for the well-known formulas 
for the transition probabilities between odd oscillator states via Legendre functions; the latter can be found in [7|,|8|, for 
example. And this really takes place, as there exist some relations between Legendre functions and Jacobi polynomials 
with specific parameter values; see [y] for details. 
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The last point to discuss is the generating functions for the transition probabilities of the singular oscillator. For 
the excitation of the ground state (the most interesting case in practice) it is easy to obtain: 

1-P X 



Go(z) = J2 Vo n z n = ( — £J , \z\ < 1/p (12) 



n=0 



Gi(*) = $>!„** = [— (- 2JP (T—) +Z ' |Z|<1/P (13) 



n=0 

while for the transitions from the first excited state (to = 1) 

The expressions (p~2]) and (fT3|) can be expanded analytically for all the complex plane. 

The formulas (|12|) and (|13J) are quite convenient to compute various operator mean-values over the transition 
probability distributions. In particular, for the adiabatic invariant I = (H) /2u> one obtains 

1+ (n — i) 1 + P / \ 

7 L = - ~ = ^ Ji , to = 0,1. 14 

7_ m — j 1 — p 

One can show [6J] that this expression remains valid for transitions from an arbitrary level. 

To summarize, we note that in our paper the transition probabilities of the singular oscillator have been calculated 
for an arbitrary frequency u)(t). To solve the problem, the modified FDM was applied and the representation theory 
for non-compact su(l, 1) algebra has been used. The final result for the transition amplitudes has been presented in 
a self-contained form, which is rather convenient for further applications. For the transitions from low-lying states 
the expressions for the generating functions have been derived and the adiabatic invariant variation at the oscillator 
evolution has been calculated. 
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